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Backscatter  Enhancement  in  Scattering 

from  Rough  Surfaces 


1.  INTRODUCTION 


Stealth  technology  has  advanced  to  the  point  where  typical  radar  target  cross  sections  are 
so  small  that  the  design  of  future  radar  systems  will  require  very  high  sensitivity.  Ihe  small 
target  cross  section  combined  with  clutter  variability  In  space  and  time  will  require  radar 
receivers  with  dynamic  ranges  of  90  dB  (over  a  1  percent  bandwidth).  Hence,  there  Is  now  a 
need  for  accurate  models  to  predict  radar  clutter  power  levels  In  various  environmental 
situations. 

As  a  result  of  recent  experimental  data  at  optical  frequencies,1  there  Is  a  question  whether 
the  traditional  clutter  prediction  models  are  accurate  under  some  conditions.  This  report 
shows  that  a  one-dimensional  physical  optics  Integral  representation  confirms  the  new  data. 
The  model  Is  then  applied  at  microwave  frequencies.  The  results  show  that  the  conventional 
model  for  mean  clutter  cross  section  can  underestimate  the  value  by  as  much  as  10  dB  for  very 
rough  surfaces. 


(Received  for  Publication  16  June  1989) 

1  Mendez,  E.R.  and  O’Donnell.  K.A.  (1987)  Observation  of  depolarization  and  backscatterlng 
enhancement  In  light  scattering  from  Gaussian  random  surfaces.  Optics  Communications, 
61(2):91-95. 
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1.1  Discussion  of  the  Experiment 


In  January  1987.  Mendez  and  O'Donnell1  first  reported  experimental  observations  of 
enhanced  backs  cattering  in  the  antispecular  direction  when  light  was  scattered  from  rough 
metallic  surfaces.  This  Initial  work  was  further  corroborated  in  a  paper2  which  reported  on 
light  scattering  from  metallic  surfaces  under  a  variety  of  rough  surface  conditions.  O’Donnell 
and  Mendez2  found  that  when  the  surface  slopes  are  small,  the  diffuse  scattering  often  agrees 
with  the  physical  optics  (PO)  theory  of  Beckmann  and  Spizzi chino3  (conventional  model)  when 
the  angle  of  incidence  is  not  too  large.  However,  when  the  surface  slopes  became  larger,  the 
theory  of  Reference  3  departs  from  the  experimental  data.  Mendez  and  O'Donnell  speculate 
that  multiple  scattering  is  a  factor. 


1.2  Summary  of  Analytical  Efforts 


There  are  a  large  number  of  theoretical  approaches4-  5-  6-  7-  8-  9- 10  to  scattering  from  rough 
surfaces.  Each  approach,  even  the  more  rigorous  theories  which  Include  multiple 
scattering,6-7  9- 10  are  based  upon  certain  restrictive  assumptions.  To  date,  only  References  8, 
9,  and  10  appear  to  predict  backscatter  enhancement  in  the  antispecular  direction  and  as 
Flood11  points  out,  of  these  only  Bahar's8  approach  shows  any  agreement  with  the 


2  O’Donnell.  K.A.  and  Mendez.  E.R.  (1987)  Experimental  study  of  scattering  from 
characterized  random  surfaces,  Journal  Optical  Society  of  America,  4:1 194-1205. 

3  Beckmann,  P.  and  Spizzi  chino,  A.  (1963)  The  Scattering  of  Electromagnetic  Waves  From 
Rough  Surfaces,  Pergamon  Press,  New  York. 

4  Desanto.  J.A.  and  Brown,  G.S.  (1986)  Analytical  techniques  for  multiple  scattering  from 
rough  surfaces.  E.  Wolf,  Progress  In  Optics,  XXIII,  Elsevier  Science  Publishers. 

5  Ruck.  G.T.,  Barrick,  D.E.,  Stuart,  W.D.  and  Krickbaum,  C.T.  (1970)  Radar  Cross  Section 
Handbook,  2,  Plenum  Press. 

6  Brown,  G.S.  (1984)  Application  of  the  integral  equation  method  of  smoothing  to  random 
surface  scattering,  IEEE  Trans,  on  Antennas  and  Propagation,  AP-32(12):  1308- 1312. 

7  Nieto- Vesperimas.  M.  and  Garcia.  N.  (1981)  A  detailed  study  of  the  scattering  of  scalar 
waves  from  random  rough  surfaces,  Opt  Acta,  28:1651-1672. 

8  Bahar,  E.  (1987)  Review  of  the  full  wave  solutions  for  rough  surface  scattering  and 
depolarization:  Comparisons  with  geometric  and  physical  optics,  perturbation  and  two-scale 
hybrid  solutions,  J.  Geophys.  Res.,  92(C5):5209-5224. 

9  Celll,  V.,  Maraduden,  A.A.,  Marvin,  A.M.  and  McGum,  A.R.  (1985)  Some  aspects  of  light 
scattering  from  a  randomly  rough  metal  surface,  J.  Opt  Soc.  America,  A2:2225. 

10  McGum,  A.R.,  Maraduden,  A  .A.  and  CeDi,  V.  (1985)  Localization  effects  in  the  scattering  of 
light  from  a  randomly  rough  grating,  Phys.  Rev.,  B3 1:4866. 

11  Flood,  W.  (1987)  Wave  Propagation:  Although  Electromagnetics  is  a  Mature  Science,  New 
Problems  are  Challenging  the  Theoreticians,  R  &  D  Roundup,  Microwaves  and  RF,  65. 
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experimental  data  of  Mender  and  O'Donnell.  This  Is  significant  in  terms  of  the  original 
experiment  since  although  multiple  levels  ■■■  roughness  appear  In  the  model.  Bahar  does  not 
include  multiple  scattering  in  hi-?  formalism. 

Thorsos  and  ishima;  u1’  have  examined  the  accuracy  of  Bahar’s®  full-wave  solution  using  a 
one-dimenstonal  sui  face  with  a  Gaussian  roughness  spectrum.  The  normalized  bistatic 
scattering  cross  section  rj°  is  c  iV-d, red  by  solving  an  iniegral  equation  for  the  scattered  field 
and  then  using  a  Monte  Carlo  technique  to  ensern.  Me  average  over  the  surface  heights.  Based 
upon  comparisons  with  this  rigorous  mathematical  simulation,  the  full- wave  solution  was 
found  to  agree  closely  with  D -  KJrehhofif  solution  '■),  except  for  small  grazing  Incident  and 
scattering  angles.  However,  Thorsos  and  fahl T«i'  ~  found  that  th?  full  ware  solution  does  not 
reduce  to  the  first  order  perturbation  solid  ion  (Ruck  st  al}3  as  the  rma  surface  height  and  slope 
become  small.  Thorsos  and  Ishimaru12  conclude  that  the  full  wave  method  is  not  an 
improvement  over  the  Klrchhoff  approx- mation  It  t  '.he  case  where  the  incident  and  scattering 
angles  are  not  close  to  90-  (grazing). 

Brown13  has  derived  an  expression  foi  'he  i.«v  ludlzd  cross  section  of  a  rough  surface  <j° 
starting  with  the  expression  for  the  far  field  approjdraatlon  for  the  scattered  field  and  making 
the  Klrchhoff  approximation  foi  the  surface  cun  eat  density  Brown  assumes  Gaussian 
distributed  surface  heights  ard  a  Gaussian  cm relation  function,  and  retains  terms  in  the 
expression  for  the  scattered  field  that  depend  upon  the  surface  slopes.  He  shows  that  for  small. 
Intermediate,  and  large  values  of  kT.  there  is  enhancement  in  the  backscatter  (anti specular) 
direction,  except  for  «.  =  0°  Here.  6t  ~  angle  of  incidence  with  respect  to  normal,  k  =  2 n/k, 

X  =  era  wavelength  and  T  =  surface  correlation  length. 

In  this  report,  a  model  for  o°  Is  used  that  is  described  in  Papa  and  Woodworth. 14  It  is 
assumed  that  the  surface  current  density  is  given  by  the  Klrchhoff  approximation.  Several 
conditions  are  placed  on  the  surface.  First,  the  surface  heights  z  are  a  random  stationary 
stochastic  process  and  are  a  function  o i  only  one  coordinate  x  (one-dimensional  roughness). 
The  heights  are  assumed  to  be  Gaussian  with  a  Gaussian  correlation.  Arbitrary  surface  slopes 
are  allowed  so  the  surface  heights  and  slopes  are  not  deeorrelated.  Under  these  conditions  the 
expression  for  a"  reduces  to  a  triple  integral,  which  may  be  evaluated  numerically.  For 
intermediate  or  larg  rtns  surface  slopes,  this  exact  physical  optics  (PO)  model  predicts 
enhanced  backseat  tc- . 


12  Thorsos,  E.I.  and  Ishimaru,  A.  (1988)  An  examination  of  the  full-wave  method  for  rough 
surface  scattering,  URSJ  Radio  Science  Session  B/F  )  dan.  1988,  Unlv.  of  Colorado,  Boulder, 
CO. 

13  Brown,  G  S  (1988)  Enhanced  Backscatter  f.ora  the  Klrchhoff  Approximation,  Private 
Communication. 

14  Papa.  R.J.  and  Woodworth,  M.B.  (1988)  T7ie  numerical  evaluation  of  a  physical  optics 
normalized  cross  sect  lot  i  for  a  rough  surface.  RADC-TR-87-280,  ADA198919! 
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2.  a°  MODELS  USED  IN  THIS  REPORT 


A  number  of  different  models  that  Incorporate  different  assumptions  are  compared  In  this 
report  In  terms  of  their  ability  to  predict  enhanced  backscatter.  This  section  will  describe 
these  models  and  their  respective  differences. 


2.1  Physical  Optics  Triple  Integral  Solution 

In  the  report  by  Papa  and  Woodworth. 14  a  general  expression  for  the  normalized  cross 
section  o°  is  derived  for  a  rough  surface  whose  random  heights,  z.  have  only  a  one- 
dlmenslonal  variation  z  *  z(x).  The  surface  current  density  is  assumed  to  be  given  by  the 
Kirchhoflf  approximation.  This  physical  optics  (PO)  result  Is  valid  only  under  the  assumption 
T  »  X  (Papa  and  Lennon).15  The  scattering  Is  assumed  to  take  place  in  the  plane  of  Incidence 
(azimuthal  angle,  <>  =  0°).  the  Incident  field  Is  horizontally  polarized  and  there  is  no 
shadowing. 

The  general  PO  expression  for  o°  for  a  one-dimensional  rough  surface  Involves  a  six-fold 
Integral  over  the  variables  xp  x2.  zp  Z2.  pp  and  p2.  where  xi  and  x2  are  two  points  on  the 
surface.  zt  and  Z2  are  the  heights,  and  pt  and  pj  are  the  slopes  at  the  two  points.  By  assuming 
that  the  random  rough  surface  Is  a  stationary  stochastic  process,  the  six-fold  Integral  reduces 
to  a  four-fold  Integral  over  the  variables  %,  x,  p{  and  p2,  where  £  ■  zj  -  Z2  and  t  ■  xj  -  x2.  If  the 
trivariate  distribution  function  for  the  heights  z  and  slopes  pi  and  p2  is  assumed  to  be 
Gaussian,  the  Integration  over  %  can  be  accomplished  analytically.  The  remaining  expression 
for  o°  then  reduces  to  a  triple  Integral  (TI).  which  may  be  evaluated  numerically  using  a 
modified  Gaussian  quadrature  technique.14 

The  expression  for  o°  was  written  as: 


°0“2fJt2jJ  dTCO*IVlW(T)  (1) 


where 


vx*  k(sln  0,  -  sin  0J 
W(t)  -  G(t)  -  H 


15  Papa.  R.J.  and  Lennon.  J.F.  (1988)  Conditions  for  the  validity  of  physical  optics  In  rough 
surface  scattering,  IEEE  Trans,  on  Antennas  and  Propagation  36(5):647-650. 
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and 


G(t)  =  2 


dp2F(P2) 


dp,  F  (p, )  |*  d4  cos(vz4)  p3  tHj.  p2,  4.  t) . 

Jo 


(2) 


Here,  p3  =  til  variate  distribution  function  In  4.  Pi  and  p2  and  v2  =  -k  (cos  0j  +  cos  0J. 

Also,  F(p)  Is  a  complicated  function  of  the  slopes  p,  the  angle  of  Incidence  0i,  the  scattering 
angle  08  and  the  Fresnel  reflection  coefficient,  and  is  given  in  Reference  14,  as  Is  the  function 
H.  There,  it  is  shown  that 


I 


I  *  I  d4  c°s(vz4)  p3  *e“c  cos(vzB) 


e-v^/4A 


(3) 


where 


M. . 

A  = - LL 


(2|R|) 
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\ 
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( 

a 

(  \ 

1 

2  2 

lMi  i  J 

Pi+Pa  +  2Fp1p2 
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-exp 


/ 

-X 

l-2t2 

V  1  V 

<J  =  standard  deviation  in  surface  heights  and  My  is  the  co-factor  of  the  covariance 
matrix  Rg.  The  function  W(x)  in  Eq.  (1)  has  singularities  in  x,  but  they  can  be  shown  to  be 
lntegrable.  This  result  forms  the  basic  expression  to  be  used  in  the  comparisons  of  this  report. 
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2.2  Small  Slope  Approximations 


Beckmann  and  Spizzl chino3  show  that  If  the  rma  surface  slopes  (V2  a/T)  are  small,  the  PO 
expression  for  a°  for  a  one-dimensional  rough  surface  reduces  to  a  single  integral  over  t: 


2krV„)  r 


4n 


1 


dt  cos!  v 


xTl 


2 

X2-Xl 


(4) 


where 


Xj,  -exp 


-oV(l-0 


;»exp 


v  A 


X. 


2  2 
V  O 

- S _ 


and  p.p  Is  the  slope  of  a  facet  that  gives  rise  to  scattering  In  the  specular  direction.  This 
simplification  is  called  physical  optics,  Ftp,p).  There  are  other,  alternative,  simplified  forms 
for  o°,  which  may  be  obtained  from  the  triple  Integral  given  by  Eqs.  (1).  (2).  and  (3).  Another 
simplification  occurs  for  a/T  <  1;  In  the  Gaussian  trl variate  distribution  function  In  Eq.  (2). 
the  heights  decorrelate  from  the  slopes  so  that 


G(*)  — >2  J*  *>p1F(|i1)ier) 


(5) 


with  I(t)  given  by  Eq.  (3).  Examination  of  this  result  shows  that  now  one  can  separately 
perform  the  remaining  Integrations  over  pi,  p2  and  t.  This  simplification  into  three  separate 
integrations  Is  called  physical  optics  with  the  F  Integral  (PO,  F  INT). 

Another  simplification  occurs  for  small  rms  slopes  (a/T  <  1);  one  argues  that  most  of  the 
contributions  to  the  pi  and  pa  integrations  In  Eq.  (5)  come  at  pi-  0  and  pa«  0  ,  so  that  then 
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G(t)  — >2  (0)  I  (x)  . 


(6) 


This  simplification  Is  called  PO,  F(0). 


2.3  Geometrical  Optica  Solution 

Let  £  -  -v^o  =  j  cos  0t  +  cos  0^, 


where  ]£  is  the  Rayleigh  parameter.  For  the  geometrical  optics  (GO)  limit.  X  -»  0  and  £  >  1  so 
that  Eq.  (4)  may  be  evaluated  for  a°  using  the  method  of  stationary  phase3 : 


o 


o 


exp 


vV^ 

43? 


(7) 


2.4  RMS  81ope  Model 

In  order  to  make  the  exact  PO.  triple  Integral  model  more  flexible  and  perhaps  improve  our 
ability  to  match  the  theoretical  model  to  experimental  data,  we  may  introduce  a  two 
correlation  lengths  surface  roughness  scattering  model,  as  follows.  In  Eqs.  (1)  through  (5),  let 
o/T  -*  o/T^hort-  where  V2  a/ T«hort  represents  the  rms  slopes  of  the  scattering  facets.  This 
implies  that,  everywhere. 


Z 

I 


fo 


2  K  (cos  0(  +  COS  ©B)| 


} 


However,  In  the  expressions  containing  the  correlation  function 
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let  T  ->  Tiong,  where  Tiong  represents  the  distance  one  must  travel  from  one  point  on  the 
surface  to  a  second  point  for  statistical  decorrelation  to  occur.  Tiong  18  the  distance  on  the 
surface  one  must  travel  from  an  initial  point  for  the  correlation  function  to  drop  to  1/e  of  its 
maximum  value.  T*hort  is  representative  of  the  rms  slopes  of  the  surface  facets  (V2  o/T*hort)- 


3.  THE  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR  METALLIC  ROUGH 
SURFACES  AT  OPTICAL  FREQUENCIES 


In  this  section,  an  attempt  will  be  made  to  compare  the  one-dimensional,  exact  triple 
Integral  solution  to  some  of  the  data  taken  at  optical  frequencies  by  O'Donnell  and  Mendez.1-2 
They  made  a  number  of  measurements  of  the  diffuse  power  scattered  from  rough  metallic 
surfaces  at  optical  frequencies  for  various  angles  of  incidence. 

Note,  after  the  figures  in  Sections  3  and  4  were  completed,  we  discovered  an  error  in  the 
factor  outside  the  integral.  It  should  be  k/(4x3)  rather  than  the  original  itk.  The  shapes  of  the 
calculated  curves  do  not  change,  but  the  magnitudes  should  be  shifted  by  a  factor  of  1/(4k9)  or 
-20.9  dB. 

In  the  experiments  of  Mendez  and  O'Donnell,1  they  scattered  a  laser  beam  (A.  *  0.633  pm) 
from  a  rough  aluminum  surface.  They  determined  the  standard  deviation  in  surface  height  o 
to  be  1.0  to  2.0  pm.  and  the  surface  correlation  function  to  be  Gaussian  with  a  1/e  width  T  -  1.8 
pm.  They  observed  that  there  is  no  backscatter  enhancement  in  the  antlspecular  (6.  «  -©J 
direction,  when  0j  =  20°.  They  state  that  the  enhanced  backscattering  can  be  easily  seen  for 
angles  of  incidence  0j  >  60°. 

In  Figure  1.  the  one-dimensional  PO  models  presented  in  Section  2  (o°  models)  are  plotted 
vs  the  scattering  angle  6a  for  an  aluminum  surface  with  e  ■  -25.3  +  J  15. 1 .  This  value  of  the 
complex  dielectric  constant  for  metals  at  optical  frequencies  is  based  upon  the  free  electron 
model  (see  Bom  and  Wolf)16  of  a  metal.  The  surface  correlation  length  is  taken  to  be  T  *  1.8 
pm.  In  the  first  sub-figure  of  Figure  1.  ■  20°,  o»1.5  pm,  and  there  is  no  backscatter 

enhancement.  This  Is  in  agreement  with  Mendez  and  O'Donnell's  data.  In  the  second  sub¬ 
figure  of  Figure  1,0!  =  20°,  o  -  2  pm,  and  the  exact  PO  triple  integral  solution  shows  no 
backscatter  enhancement.  In  the  third  sub-figure  of  Figure  1.  0j  -  20°  and  the  rms  surface 
slopes  have  been  greatly  increased  so  o/T  -  2.78;  however,  there  is  a  slight  backscatter 
enhancement  at  0.  -  -90°.  In  the  fourth  sub-figure  of  Figure  1 ,  ©t  has  been  increased  to  50°  and 
o/T  -l.ll.  It  may  be  readily  observed  that  there  is  backscatter  enhancement  at  0,  -  -90°. 

This  also  is  in  agreement  with  Mendez  and  O'Donnell's  data.  It  should  be  noted  that  the  PO, 
Ffpap)  approximation  and  the  GO  approximation  also  predict  backscatter  enhancement  in  the 
third  and  fourth  sub-figures  of  Figure  1,  but  these  results  are  not  reliable;  both  approximations 
are  beyond  their  range  of  validity  since  they  were  derived  under  the  assumptions  of  small  rms 
surface  slopes,  (o/T  <  1). 


16  Bom.  M.  and  Wolf,  E.  (1959)  Principles  of  Optics,  Pergamon  Press. 
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Note,  quite  often  these  two  approximations  predict  the  same  RCS  a°  values,  hence  the  lines 
are  plotted  on  top  of  one  another,  giving  the  appearance  of  a  long  dashed  bne. 

In  Figure  2,  the  angle  of  Incidence  0i  Is  equal  to  70°  and  ols  1  pm  with  T  «  1.8  pm  for  a 
rough  aluminum  surface  Irradiated  with  laser  light  (X  a  0.633  pm).  These  are  the  parameters 
given  by  Mendez  and  O'Donnell.1  It  may  be  readily  discerned  that  the  triple  Integral  solution 
predicts  enhanced  backscatter,  and  also  by  the  PO,  F(pgp)  model  and  the  GO  model.  For 
0,  =  70°.  Mendez  and  O'Donnell  reported  observation  of  enhanced  backscatter  for  aluminum. 
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et  70.00 
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Figure  2.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  an  Aluminum  Surface,  0i  =  70°,  a/T  -  0.556. 

In  Figure  3,  all  the  parameters  are  the  same  as  In  Figure  2,  except  that  the  upper  limit 
for  o  reported  by  Mendez  and  O'Donnell  Is  used  (a  =  2  pm).  Here,  again,  for  8i  =  70°,  the  triple 
Integral  solution  predicts  enhanced  backscatter,  In  accord  with  Mendez  and  O'Donnell's  data. 
Once  again,  the  PO,  F(p^p)  model  and  the  GO  model  predict  enhanced  backscatter,  even  though 
both  models  are  being  used  beyond  their  range  of  validity  ( a/T  =  1.11). 
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Figure  3.  The  Normalized  Cross  Section  a°  vs  Scattering  Angle 
for  an  Aluminum  Surface,  0j  =  70°,  o/T  =  1.11. 

In  Figure  4  the  parameters  are  the  same  as  in  Figure  3.  except  that  the  angle  of  Incidence 
0i  =  60°.  Here,  also,  there  is  enhanced  backscatter,  although  it  is  not  as  pronounced  as  In 
Figure  3  (0t  =  70°).  This  trend  Is  in  agreement  with  the  experimental  data. 
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Figure  4.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  an  Aluminum  Surface,  »  60°,  o/T  =  1.11. 

In  Figure  5,  9i  ■  70°  and  the  correlation  length  for  a  rough  aluminum  surface  has  been 
Increased  to  T  -  4.5  pm.  Again,  It  may  be  readily  discerned  that  there  Is  enhanced  backscatter. 
This  Is  also  true  In  Figure  6,  where  0j  ■  70°,  T  ■  1 .8  pm  and  the  standard  deviation  In  surface 
height  has  been  decreased  too*  0.9  pm.  Figure  7  shows  the  normalized  received  power 
scattered  from  a  gold  coated  diffuser  at  X  *  0.6333  pm  with  a  vertically  polarized  incident 
wave.  The  normalized  Intensity  (In  dB)  is  plotted  vs  the  scattering  angle  for  +  -  0°  ( in  the 
plane  of  Incidence).  O'Donnell  and  Mendez3  formed  a  histogram  (Gaussian)  of  the  surface 
height  data  measured  with  a  profilometer  the  standard  deviation  in  height  was  o  ■  2.27  pm. 
The  surface  correlation  function  was  also  measured  and  found  to  be  Gaussian  with  a 
correlation  length  T  ■  20.9  pm.  The  free  electron  theory  of  metals  Irradiated  at  optical 
frequencies16  gives  the  following  value  for  the  complex  dielectric  constant  for  gold: 
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e  «  -7.8  +  j  2.7.  It  may  be  discerned  from  Figure  7  that  O'Donnell  and  Mendez's2  experimental 
data  clearly  shows  enhanced  backscatter  in  the  andspecular  direction,  0„  *  -70°.  In  Figure  8, 
the  o°  values  are  plotted  vs  0,  for  the  triple  integral  solution  and  the  various  approximations 
for  the  conditions  of  O'Donnell  and  Mendez's2  gold  diffuser  experiment.  It  is  readily  observed 
that  none  of  the  theoretical  models  predict  enhanced  backscatter.  This  is  because  all  these 
models  are  single  scattering  theories,  and  provide  no  physical  mechanism  for  enhanced 
backscatter  at  such  low  rms  surface  slopes  (o/T  *  0.1).  When  the  rms  surface  slopes  are 
intermediate  or  large  (o/T  2  1),  then  there  are  enough  facets  oriented  so  that  specular 
scattering  from  these  facets  can  reach  a  maximum  in  the  backward  direction.  In  this  case,  the 
single  scatter  physical  optics  model  can  show  backscatter  enhancement. 
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Figure  5.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  an  Aluminum  Surface,  0j  *  70°,  o/T  ■  0.5. 
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Figure  6.  The  Normalized  Cross  Section  <j°  vs  Scattering  Angle 
for  an  Aluminum  Surface.  0i  =  70°,  a/T  =  0.5  with  Large 
Standard  Deviation  in  Surface  Height. 
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Figure  7.  The  Normalized  Intensity  vs  Scattering  Angle  for  a 
Gold  Diffuser,  Experimental  Data 
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Figure  8.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface,  0j  ■  70°,  o/T  •  0.109. 

In  Figure  9,  t  ■  -7.8  +  j  2.7  for  gold,  0j  -  70°  and  the  rms  slope  Is  Increased  from 
o/T  »  0.109  to  0.216  by  decreasing  the  correlation  length,  T  *  10.5  pm.  Here,  the  rms  surface 
slopes  are  still  too  small  to  produce  any  enhanced  backscatter.  In  Figure  10,  the  parameters 
are  the  same  as  in  Figure  9,  except  for  the  correlation  length,  which  has  been  further  reduced 
to  T  m  5.25  pm.  Now.  o/T  -  0.432  is  sufficiently  large  so  that  there  Is  a  small  amount  of 
enhanced  backscatter.  In  Figure  11.  the  parameters  are  the  same  as  In  Figure  10,  except  for  the 
correlation  length,  which  again  has  been  even  further  reduced  to  T  ■  2.0  pm.  It  may  easily  be 
discerned  that  with  o/T  ■  1.135,  there  Is  considerable  backscatter  enhancement  for  this  gold 
diffuser. 


17 


SCATTER  ANGLE  (deg) 


e,  70.00 

0.633E-06 
o  0.227E-05 

T  0. 105E-04 

o/T  0.216 

e  -7.8  2.7 


-  TRIPLE  ENT 

- PO.  FINT 

- PO.  F(H.p) 

-  PO.  F(0) 

- GO,  Ft  (Xap ) 


Figure  9.  The  Normalized  Croes  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface,  0j  ■  70s,  o/T  ■  0.216. 
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Figure  10.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface,  0|  =  70°,  a/T  =  0.432. 
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Figure  11.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface,  0,  ■  70°.  a/T  m  1.135. 

In  Figures  12,  13.  and  14.  the  two  correlation  lengths  surface  roughness  model  discussed 
in  Section  2  is  investigated.  Here,  e  -  -7.8  +  j  2.7  for  a  gold  plated  rough  surface  and 
X  m  0.633  pm.  The  correlation  length  that  governs  the  decorrelation  properties  of  the 
scattering  facets  Ti^  -  20.9  pm  in  these  three  figures,  in  accordance  with  O'Donnell  and 
Mendez's2  experimental  data.  The  correlation  length  T»hort  that  governs  the  rms  surface  slopes 
in  Figure  12  is  chosen  so  that  o/T,hort  *  0.454.  It  may  be  readily  noted  that  there  is  no 
backscatter  at  all  from  such  a  surface,  much  less  backs catter  enhancement.  In  Figure  13, 

Tshort  has  been  reduced  so  that  o/T*kart  ■  1.14.  However,  even  with  an  rms  surface  slope 
greater  than  1,  there  is  no  backscatter  enhancement.  In  Figure  14,  Tahort  has  been  further 
reduced  so  that  now  o/Tahort  ”  2.27.  Here,  one  may  observe  that  there  is  some  backscatter.  but 
not  enough  to  be  considered  an  enhancement  in  the  backscatter  direction. 
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Figure  12.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface,  0i  =  70°,  Two  Correlation  Lengths  cr/TShort  =  0.4504. 
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Figure  13.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface.  0j  =  70°,  Two  Correlation  Lengths  ff/T,hort  =  1. 14. 
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Figure  14.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Gold  Surface.  0i  =  70°.  Two  Correlation  Lengths  o/T8h0rt  =  2.27. 


4.  THE  THEORETICAL  RESULTS  FOR  DIELECTRIC  ROUGH  SURFACES  AT  RADIO 

FREQUENCIES 

For  metallic  rough  surfaces,  the  theoretical  model  for  o°  predicts  backscatter  enhancement 
for  moderate  to  large  rms  surface  slopes,  in  agreement  with  some  experimental  data.  This 
would  suggest  that  backscatter  enhancement  might  occur  for  dielectric  rough  surfaces  also.  In 
this  section,  the  normalized  scattering  cross  section  a°  will  be  studied  for  rough  dielectric 
surfaces  as  a  function  of  elevation  scattering  angle  9,  for  various  rms  surface  slopes  V2  o/T 
and  angles  of  incidence  0t.  Comparisons  will  be  made  between  the  exact  PO,  triple  integral 
solution  and  various  approximations:  the  PO,  F  INT,  FCpap),  F(0)  and  the  high  frequency  GO 
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model.  Also.  It  will  be  shown  that  for  moderate  to  high  rms  surface  slopes  (o/T  2  1),  there  Is 
backscatter  enhancement  near  the  antlspecular  direction  8.  -  - 

In  Figure  15,  o°  is  plotted  vs  elevation  scattering  angle  for  a  dielectric  surface  with  e  »  4.0, 
correlation  length  T  =  0.5m  and  an  Incident  em  wavelength  of  0.1m.  In  the  first  three  sub¬ 
figures,  the  angle  of  Incidence  0i  a  75°,  and  the  rms  surface  slopes  successively  Increase  from 
o/T  *  0.2,  to  1.0  up  to  5.00.  It  may  readily  be  observed  that  for  small  rms  slopes,  (o/T  *  0.2), 
there  Is  no  enhanced  backscatter.  There  Is  good  agreement  between  the  exact  triple  Integral 
solution  and  the  PO,  Ffyisp)  model  and  the  GO  model.  The  PO,  F  INT  model  and  the  PO,  F(0) 
model  show  some  discrepancy  when  compared  with  the  exact  triple  Integral  solution.  In  the 
second  sub-figure,  o/T  =  1 .0,  and  only  the  triple  integral  solution  shows  enhanced  backscatter. 
The  PO,  Ffyisp)  model  and  the  GO  model  show  no  enhanced  backscatter,  the  PO,  F  INT  model 
and  the  PO,  F(0)  model  show  large  discrepancies  with  the  triple  Integral  solution.  In  the  third 
sub-figure,  o/T  =  5.0,  and  all  the  approximate  models  show  discrepancies  with  the  exact  triple 
Integral  solution.  The  triple  Integral  solution  shows  enhanced  backscatter  at  0S  -  -90°.  In  the 
fourth  sub-figure  of  Figure  15,  the  angle  of  Incidence  «  89.75°  and  o/T  *  1 .0.  Here  again,  the 
triple  Integral  solution  shows  enhanced  backscatter,  with  all  the  approximate  models 
exhibiting  discrepancies  compared  with  the  triple  Integral  solution. 
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Figure  15.  The  Normalised  Cross  Section  o°  vs  Scattering  Angle 
for  a  Dielectric  Surface  with  Different  Roughness  Conditions. 
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In  Figure  16.  the  angle  of  incidence  is  reduced  to  0j «  60°  and  the  dielectric  constant  e 
remains  4.0.  Here,  the  rms  surface  slope  has  an  intermediate  value  a/ T  •  1.0  and  the  em 
wavelength  remains  at  0.1m.  It  may  readily  be  discerned  that  the  exact  triple  integral  solution 
exhibits  the  most  enhanced  backscatter  and  the  approximate  PO  models  show  discrepancies 
compared  with  the  triple  Integral  solution. 


SCATTER  ANGLE  (deg) 


60.00 


0.100 


0.500 


0.500 


1.00 


4.00 


TRIPLE  INT 
PO.F(n.p) 
PO.  F(0) 

PO,  F  INT 


Figure  16.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Dielectric  Surface,  0j  -  60°. 


In  Figure  17,  the  normalized  cross  section  o°  of  a  rough  dielectric  surface  is  plotted  vs 
scattering  angle  0.  for  0j  *  75°  using  a  two-dimensional  PO  scattering  model,  with  a  single 
Integral  representation  for  o°  analogous  to  Eq.  (4).  It  may  be  noted  that  there  is  backscatter 
enhancement,  but  it  is  not  nearly  as  pronounced  as  the  triple  integral  solution  (in  the  one- 
dimensional  model)  shown  in  the  second  sub-figure  of  Figure  1. 
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Figure  17.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Dielectric  Surface,  0j  «*  75°. 


In  Figures  18  through  22,  the  dielectric  constant  has  been  changed  to  e  -  24.0  +  J  32.  which 
is  representative  of  wet  sandy  loam.  The  em  wavelength  is  X  ■  0.03m  (X-band).  In  Figures  18 
through  21.  the  angle  of  incidence  8|  ■  40°.  and  the  rms  surface  slopes  became  progressively 
larger.  It  may  be  noted  from  Figures  18  and  19  that  there  is  no  backscatter  enhancement  for 
o/T  <  1.0.  Figures  20  and  21  show  that  there  Is  backscatter  enhancement  for  o/T  >  1.0,  but 
that  it  occurs  near  0,  *  -90°. 
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Figure  2 1 .  The  Normalized  Cross  Section  cr°  vs  Scattering  Angle 
for  a  Lossy  Dielectric  Surface,  0i  =  40°.  o/T  =  1.67. 
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Figure  22.  The  Normalized  Cross  Section  o°  vs  Scattering  Angle 
for  a  Lossy  Dielectric  Surface,  0!  »  75°,  o/T  ■  1.25. 
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In  Figure  22.  the  angle  of  Incidence  has  been  Increased  to  0i  ■  75°  for  the  dielectric  of 
e  ■  24  +  J  32.  With  a  moderate  rms  slope  (o/T  ■  1 .25),  It  may  be  noted  that  the  triple  integral 
solution  and  the  PO,  F(p.p)  approximation  and  the  GO  model  all  show  very  large  enhanced 
backscatter.  This  is  analogous  to  the  behavior  for  metallic  surfaces,  which  is  to  be  expected 
since  a  large  Imaginary  part  of  the  complex  dielectric  constant  causes  the  dielectric  to  behave 
more  like  a  metallic  material.  The  two  PO  approximations,  F(0)  and  F  INT  show  very  poor 
agreement  with  the  TI  solution.  It  is  to  be  expected  that  the  PO,  F(0)  approximation  shows  veiy 
poor  agreement  with  the  exact  triple  Integral  solution  when  o/T  >  1,  since  It  was  derived  under 
assumptions  that  are  more  restrictive  than  the  Ffyt*p)  and  GO  models;  the  rms  surface  slopes 
are  not  Just  small,  o/T  must  be  approximately  zero  (o/T  -  0). 
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S.  SUMMARY  AND  CONCLUSIONS 


In  this  report,  It  has  been  demonstrated  that  backs catter  enhancement  can  occur  for  both 
rough  dielectric  surfaces  at  microwave  frequencies  and  metallic  rough  surfaces  at  optical 
frequencies.  This  Is  observed  only  when  the  angle  of  Incidence  Is  greater  than  or  equal  to  50° 
and  the  rms  slope  of  the  surface  Irregularities  Is  equal  to  or  greater  than  1 .0.  The  backscatter 
enhancement  Is  predicted  from  a  physical  optics  (PO)  model  that  is  based  on  a  rough  surface 
scattering  formalism  In  which  it  is  assumed  that  the  surface  height  z  varies  only  In  one 
transverse  dimension  x.  The  relatively  large  rms  surface  heights  Imply  that  the  multivariate 
distribution  function  describing  the  statistics  of  the  surface  heights  and  slopes  do  not 
decorrelate.  This  results  In  a  formula  for  the  normalized  cross  section  o°  for  the  rough  surface 
which  must  be  expressed  as  a  triple  integral  fn).  Simpler  PO  models  can  be  obtained  by 
assuming  the  rms  surface  slopes  are  small  (<j/T  <  1).  Still  further  simplifications  for  o°  can  be 
obtained  by  taking  the  high  frequency  geometrical  optics  (GO)  limit.  The  various  models  for  o° 
are  compared  for  various  angles  of  Incidence  0(  and  scattering  6S  as  a  function  of  surface 
roughness.  In  general,  only  the  exact  H  model  accurately  predicts  enhanced  em  wave 
backscatterlng  near  the  antlspecular  (0.  =  -  0J  direction.  Of  the  various  approximations  for  o°, 
the  single  Integral  PO  F(p^,)  model  agrees  best  with  the  exact  TI  solution  for  most  ranges  of  the 
parameters  (0t,  0»,  o  and  T).  In  the  examples  shown  here,  GO  also  gives  good  agreement. 

The  physical  mechanism  responsible  for  the  enhanced  backscatter  predicted  by  the  PO  TI 
model  is  not  multiple  scattering,  because  PO  Is  based  on  a  single  scatter  formalism.  Here,  it  is 
diffuse  backscatter  from  a  sufficient  number  of  facets  with  relatively  large  rms  surface  slopes. 
This  Is  the  reason  why  the  TI  model  gave  good  agreement  with  the  experimental  results  of 
Mendez  and  O'Donnell,  where  the  rms  slopes  for  the  aluminum  surface  were  relatively  large. 
However,  for  the  scattering  from  a  rough  gold  surface  reported  by  O'Donnell  and  Mendez2,  the 
agreement  between  the  theoretical  triple  integral  model  and  the  experimental  data  was  not  too 
good  because  the  rms  surface  slopes  were  small  (of T  <  0.5).  O'Donnell  and  Mendez2  claim  that 
in  this  case  the  backscatter  enhancement  must  be  due  to  multiple  scattering.  Therefore, 
although  multiple  scattering  may  Increase  as  the  rms  slope  increases,  a  second  mechanism 
comes  Into  play  for  moderate  and  large  slopes,  that  Is,  there  are  enough  facets  oriented  to  give 
a  specular  return  in  the  backscatter  direction. 
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The  theoretical  TI  model  will  predict  backscatter  enhancement  for  a  gold  surface  If 
o/T  >  0.8.  Bahar17  has  observed  backscatter  enhancement  at  small  rms  surface  slopes  using 
the  full  wave  solution,  so  that  multiple  scattering  Is  not  the  only  mechanism  here.  Bahar' s 
observation  of  backscatter  enhancement  was  more  pronounced  when  a  Plerson-Moskowltz18 
spectrum  was  used  to  describe  the  correlation  properties  of  the  surface,  rather  than  a  Gaussian 
spectrum.  Bahar17  claims  that  the  Plerson-Moskowltz  spectrum  better  fits  the  correlation 
properties  of  the  surfaces  measured  by  Mendez  and  O'Donnell.12 


17  Bahar,  E.  (1989),  Private  Communication. 

18  Brown,  G.  (1978)  Backscatter  from  a  Gaussian-distributed  perfectly  conducting  rough 
surface,  IEEE  Trans,  on  Antennas  and  Propagation  26(3):472-482. 
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